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Abstract 

In  order  to  provide  computationally  affordable  predictive  control  algorithms,  predictive  switching  logic 
schemes  are  considered  whereby  a  feedback-gain  is  switched-on  at  any  time  from  a  family  of  candidate 
feedback- gains  so  as  to  control  a  discrete-time  input-saturated  LTI  system  possibly  subject  to  persistent 
bounded  disturbances  of  unknown  arbitrary  magnitude.  It  is  constructively  shown  that  such  schemes  do 
exist  which  ensure,  along  with  good  tracking  performance,  global  asymptotic  and  semi-global  exponential 
stability  in  the  noiseless  case,  as  well  as  finite  loo-induced  gain  to  the  disturbance-to- state  map,  whenever 
the  structure  of  the  disturbed  plant  can  make  such  properties  conceptually  achievable. 

Key  words:  Switching  control;  Predictive  control;  Control  of  input-saturated  systems;  Anti-windup  con¬ 
trol;  Nonlinear  control. 

1  Introduction 

In  recent  years,  control  of  input-saturated  dynamical  systems  -  a  subject  of  ever-lasting  fundamental 
interest  in  control  engineering  -  has  attracted  significant  research  efforts.  These  contributions  have  been 
basically  of  a  two-fold  nature:  on  one  side  (Sussmann,  Sontag  &  Yang,  1994;  Hou,  Saberi,  Lin  &  Sannuti, 
1998)  the  major  emphasis  has  been  on  basic  issues,  e.g.,  characterization  of  state  domains  whose  elements 
can  be  asymptotically  regulated  to  zero  by  bounded  inputs;  on  the  other  side,  the  main  interest  has 
been  on  anti- windup  control  schemes  (Hanus,  Kinnaert  &  Henrotte,  1987;  Glattfelder,  Ohta,  Mosca  & 
Wieland,  2000)  whereby  suitable  corrections  to  a  pre-designed  compensator  are  generated  whenever  input 
saturations  take  place.  A  common  shortcoming  of  these  contributions  is  that,  though  at  times  robustness 
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against  model  inaccuracies  and  asymptotically  vanishing  disturbances  are  considered,  little  is  known  on 
how  to  deal  efficiently  with  persistent  disturbances  of  unknown  arbitrary  magnitude. 

A  third  source  of  contributions  to  the  subject  has  been  originated  within  model-based  predictive  control 
(MBPC)  (Mosca,  1995;  Mayne,  Rawlings,  Rao  &  Scokaert,  2000;  Maciejowski,  2002).  MBPC,  though 
inherently  tailored  to  handle  input  and/or  state-related  constraints,  has  been  so  far  hampered  in  having 
a  significant  impact  on  applications  other  than  slow-process  control,  mainly  because  of  its  high  computa¬ 
tional  load.  Moreover,  an  extension  of  MBPC  to  the  case  of  persistent  disturbances  (Scokaert  &  Mayne, 
1998)  appears  to  be  so  computationally  demanding  that  it  can  be  yet  considered  in  practice  an  unsolved 
problem. 

The  lecture  related  to  this  paper  aims  at  describing  a  computationally  affordable  solution  to  the  regu¬ 
lation  problem  of  discrete-time  input-saturated  linear  time-invariant  (LTI)  systems  subject  to  persistent 
bounded  disturbances  of  unknown  arbitrary  magnitude  .  The  solution  that  will  be  elaborated  hereafter 
enjoys  the  following  features: 

1.  It  is  realized  via  a  supervisory  switching  control  scheme  whereby  a  feedback-gain  matrix,  selected 
from  a  finite  family  of  pre-designed  candidate  feedback-gains,  is  at  any  time  switched-on  in  feedback 
to  the  plant  according  to  the  previous  gain  matrix  and  the  information,  either  complete  or  partial, 
on  the  current  plant  state; 

2.  No  disturbance  upper-bound  need  to  be  known; 

3.  The  feedback-gain  selection  is  made  in  accordance  with  a  predictive  control  philosophy,  and  each 
candidate  feedback-gain  is  tuned  on  to  a  different  control-horizon; 

4.  The  supervisory  switching  logic  is  flexible  enough  so  as  to  enable  the  designer  to  simplify  the  scheme 
by  trading  off  performance  vs.  memory  and/or  computational  complexity,  while  retaining  guaranteed 
stability  properties. 

For  an  overview  on  hybrid  systems  -  the  general  framework  which  encompasses  the  developments  of  the 
present  paper  -  the  reader  is  referred  to  (Antsaklis,  2000). 


2  Problem  Formulation 


Some  notations  used  throughout  the  paper  are  listed  hereafter  for  the  reader’s  convenience: 

m  :=  (1, 2,  •  •  •  ,  m} 

with  m  an  integer  greater  than  or  equal  to  1;  ||x||/j  is  a  shorthand  notation  for  \x'  P(h)x\1/2  where  the 
prime  denotes  transpose  and  P(h )  is  as  in  (17);  finally,  ||x||  :=  [ x'x ]1//2  denotes  the  usual  (un-weighted) 
Euclidean  norm. 

The  paper  deals  with  the  regulation  problem  of  a  discrete-time  input-saturated  LTI  plant 

f  x(t  +  1)  =  $x(t)  +  Ga(u(t))  +  £(f)  ,  s 

\  x(t)  =  x(t)  +  x{t) 


where:  t  E  Z+  :=  {0, 1, . . x  E  Mn  is  the  plant  state;  u  =  [ui, 
denotes  transpose; 


u  ,  Ui^y-u^Ui 
n(u )  ,  elsewhere 


•  ■  •,  um]7  E  Mm  the  control  input;  the  prime 


Vi  E  m 


(2) 
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Ui,Ui  >  0;  and  n(u)  :  Mm  — >  Mm,  an  arbitrary  bounded  nonlinear  function  of  u.  In  (1),  both  £  and  x  are 
bounded  and  possibly  persistent  disturbances  of  unknown  arbitrary  magnitude.  The  problem  is  to  find, 
based  on  x,  a  partial  state-information  vector,  feedback  controls 

u  =  f{x )  (3) 

so  as  to  ensure,  under  suitable  conditions,  stability  in  the  noiseless  case  as  well  as  finite  loo-induced  gain 
of  the  disturbance-to-state  map  from  £  and  x  to  x  embedded  in  (l)-(3).  The  adopted  approach  consists 
of  selecting  a  discrete  family  T  =  {T)l}^=1  of  linear  state-feedback  gains  F] l  and  a  switching  logic 

h(t)  =  £(x{t),  h(t  —  1))  (4) 

in  such  a  way  that  the  regulated  plant 

x{t  +  1)  =  4 >x(t)  +  Ga(Fh(t)x(t))  +  £(i)  (5) 

have  the  stated  stability  properties. 

The  paper  is  organized  as  follows.  Sect.  3  describes  the  specific  type  of  feedback-gain  matrices  that 
are  adopted  to  realize  possible  control  actions.  Sect.  4  shows  that,  thanks  to  the  type  of  candidate 
feedback-gain  matrices,  stability  under  arbitrary  switching  control  can  be  ensured  to  an  LTI  system  with 
no  input  saturations  by  only  imposing  a  simple  but  essential  admissibility  condition  on  the  switching 
sequences.  Sect.  5  exploits  the  stability  results  in  Sect.  4  so  as  to  extend  them,  via  the  adoption 
of  appropriate  switching  logic  supervisors,  to  LTI  systems  subject  to  input  saturations  and  persistent 
bounded  disturbances  of  unknown  arbitrary  magnitude.  Sect.  6  reports  an  examples,  while  Sect.  7  ends 
the  paper  with  conclusive  remarks. 

3  Candidate  Receding  Horizon  Feedback- Gains 

Consider  temporarily  the  noiseless  linear  variant  of  (1) 

x(t  +  1)  =  4>x(t)  +  Gu(t)  (6) 

Assume 

E  =  (4>,  G)  reachable  (7) 

Note  that  (7)  entails  no  loss  of  generality  in  that,  if  E  is  stabilizable  but  not  reachable,  all  subsequent 
developments  apply  to  the  reachable  subsystem  of  the  plant  obtained  via  a  Gilbert-Kalman  reachability 
decomposition.  Because  of  properties  that  are  motivated  next  in  some  detail,  the  candidate  feedback-gains 
are  chosen  as  follows 

Fh  :=  1G'{<S>h-1)lg-h1§h  (8) 

where:  h  is  a  positive  integer,  h  >  u,  u  the  reachability  index  of  E;  4^  =  4/^  >  0;  and  Qh  is  the  h- order 
reachability  Gramian 

h 

gh  ■=  &k-1G*-1G'($k-1y  (9) 

k= 1 

An  explanation  on  where  (8)  stems  from  is  in  order.  To  this  end,  we  consider  the  following  control  problem 

Vh. 
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Zero  terminal-state  minimum  energy  control  problem  Vh-  Let  x  be  the  state  of  (6)  at  time  0, 
and  Uh{x )  the  set  of  all  controls  oj  of  length  h.  u  =  [V(0)  •  •  •  u'(h  —  1)];  which  drive  the  system  state  x  to 
the  zero-state  Ox  in  h  time-steps 

Uh{x)  :=  {w  G  (Rm)h  :  x(h)  =  Ox}  (10) 

where  x(h)  =  <&hx  +  Ylk=o  ^h~1~leGu(k).  Note  that  Uh{x )  /  0,  \/x  G  Mn,  if  h  >  v.  Vh  consists  of  finding 
Uh(x),  the  element  'mUh(x)  of  minimum  energy 


h- 1 

j2u'muU(k)  =  uj'^uoj  (ii) 

k= 0 

where  ^/u  =  ^'u  >  0  and  :=  Diag{'Pix, . . . ,  ’Lli}-(/i-times).  For  h  >  v,  such  an  element  Uh(x)  is  as 
follows 


Uh(x) 


Th 


[u'h{0\x), . . .  ,u'h(h  -  l\x)]' 

[F'h{v)---F'h(h-i)}'x 

Thx 


where  IZh  is  the  h- order  reachability  matrix 


P-h  '■= 


$h-1G|...  |$G|G 


The  integer  h  is  referred  to  as  the  control  horizon  associated  to  7 


(12) 

(13) 


(14) 


Note  that  in  (8)  is  given  in  terms  of  Fh  as  follows 

Ph  \jm  0mxm(/i_l)]  J~h  (lb) 

Hence,  Ff t  is  recognized  to  be  the  feedback-gain  matrix  of  the  receding  horizon  regulator  (Kwon  &  Pearson, 
1975)  associated  to  problem  TV 
From  (12)  and  (13)  it  follows  that 

v!h(x)$uuh(x)  =  \\x\\l  :=  x'P(h)x  (16) 

where 

p{h)  =  (^hyg^h 

=  P'{h)  >  0  (17) 

It  can  be  seen  (Mosca  1995)  that,  if  P(v)  is  as  in  (17)  for  h  =  n,  then  for  h  >  u,  P(h  +  1)  satisfies  the 
Riccati  difference  equation 

P{h  +  1)  =  &P(h)§  -  &P(h)G  (tfu  -I-  G'P(h)G) _1  G'P(/i)$  (18) 

and 

Fh+1  =  -[ttu  +  G'P{h)G]-lG'P{h)<5>  (19) 
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In  addition, 

P(h  +  1)  <  P[h) 

From  (17)  and  (20),  it  follows  that  the  following  limit  exists 

lim  P(h)  =:  P( oo)  >  0 

h—>  oo 


(20) 


(21) 


Further,  if  (6)  is  ANCBI  (asymptotically  null-controllable  with  bounded  input),  viz.  (Sussmann  et  al, 
1994),  (6)  is  stabilizable  and  has  no  exponentially  unstable  eigenvalue, 


P(  oo)  =  0n 


(22) 


with  rate  of  convergence  faster  than  or  equal  to  1/h. 

Let 

Mh(x)  :=  max{  [Uh(k\x)]i  .  y.  _  — ;  k  +  1  £  h;  i  £  rh\  (23) 

Vi 

where  [u]i  denotes  the  i-th  component  of  the  vector  u.  Note  that  the  whole  sequence  Uh(x )  stays  within 
the  linear  range  of  cr(-)  if  and  only  if  Mh(x)  <  1.  It  follows  from  (16)  and  (22)  that  for  an  ANCBI  system 
for  any  x  £  Mn  it  is  possible  to  find  a  large  enough  horizon  h  so  as  to  satisfy  M^x)  <  1. 


The  considerations  in  the  remark  that  follows  play  a  significant  role  in  the  proofs  of  Theorem  3  and 
Theorem  4. 


Remark  1  -  Consider  the  orthogonal  decomposition 

Mn  =  ft(($n)')  ©  A f(<S>n)  (24) 

where  A 7(<3?n)  =  Af(&h),  Wi  >  re,  and  TZ(-)  and  A f{-)  denote  range-space  and,  respectively,  null-space. 
Accordingly,  let  x  =  xr  +  xL ,  with  xr  £  lZ(($>ny)  and  x'rx L  =  0.  Notice  that  ||x||/t  =  ||xr||ft. 

If  =  S,  S  a  stability  matrix,  as  h  — >  oo  one  has 

UXs\\xr\\  <  Mh(x)  <  U  A|||xr||  (25) 

where  U,  U  £  (0,  oo),  and  A s  (As)  denotes  the  maximum  (the  nonzero  minimum)  among  the  moduli  of 
the  5-eigenvalues. 

Let  $  =  Q,  where  Q  has  all  its  eigenvalues  on  the  unit  circle.  From  (13)  and  the  fact  that  Q  is  similar  to 
a  upper  triangular  matrix  of  Jordan  blocks  associated  to  eigenvalues  of  unitary  modulus,  it  follows  that, 
as  h  — *  oo, 

Mh(x)  =  V(x)[l  +  0(h-1)}h-pW\\x\\  <Mh-l\\x\\  (26) 

where  0(h~l)  stands  for  a  quantity  of  the  order  of  Z©1  or  that  vanishes  at  a  faster  rate,  while  both  V(x) 
and  p(x)  only  depend  on  x/||x||.  Further,  V(x)  and  M  are  positive  reals,  and  p(x)  an  integer  which  takes 
values  in  /j  ,  where  p  denotes  the  maximum  among  the  dimensions  of  the  Jordan-blocks  associated  to  the 
Q-eigenvalues.  ■ 


In  the  remaining  part  of  this  paper,  particular  attention  will  be  devoted  to  systems  algebraically  equivalent 
to  the  ones  with  transition-matrices  having  the  following  block-diagonal  structure 


x(t  +  1) 


S  0 
0  Q 


s(t) 

Q(t) 


+ 


Gq 


cr(u(t))  + 


Ut)  ' 
£?(*)  . 


(27) 
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where:  x  =  [s'  </]';  dims  =  dimS ;  £.,(■)  and  £g(-)  are  bounded  disturbances;  S  a  stability  matrix;  and  Q 
has  all  its  eigenvalues  of  unit  modulus  with  arbitrary  geometric  multiplicities.  The  relevance  of  (27)  is 
that  a  generic  ANCBI  system  is  algebraically  equivalent  to  a  system  with  a  block-diagonal  structure  as 
(27)  with  ( Q,Gq )  reachable.  A  class  of  ANCBI  systems,  sometimes  considered  in  the  literature  (Hou  et 
al.,  1998;  Saberi  et  al. ,  2000),  is  the  one  wherein  the  stable  modes  are  affected  by  the  neutrally  stable 
ones,  while  the  latter  are  only  affected  by  the  controls,  and  the  disturbances  only  enter  the  stable  ones 


z(t  +  1)  = 

A  A 

0  A 

zs(t ) 

_  zq(t) 

+ 

'  B  ' 
B 

1 - 1 

o 

3 

where:  z  =  [z's  zq]f\  dimzs  =  dimA ;  us(t)  a  disturbance;  A  a  stability  matrix;  and  A  has  all  its  eigenvalues 
of  unit  modulus  with  arbitrary  geometric  multiplicities.  Systems  of  the  form  (28)  can  be  shown  to  be 
algebraically  equivalent  to  systems  as  in  (27)  with  dimS  =  dimA  and  £q(t)  =  0:  a  particular  instance 
considered  in  Theorem  3. 


Remark  2  -  Assume  that  the  linear  (unsaturated)  variant  of  (27)  be  reachable  by  the  input  o{u)  =  u. 
Because  (12)  is  linear  in  x, 


and,  by  (27), 


uh(x)  =  ush(s)  +  uqh(q) 


(29) 


ush(k\s)  :=  -^-1G'{^h-k-1yg^1Shs  (30) 

where  k  +  1  £  h  ,  denotes  the  state-transition  matrix  of  (27),  and  similar  definition  for  uqh(q). 

Let  ||g||  =  0.  In  such  a  case,  lower  and  upper-bounds  similar  to  the  ones  in  (25)  can  be  established,  viz. , 

U  A§||sr||  <  Mh([s'  O']')  <  U  As||sr||  (31) 


where  sr  is  the  orthogonal  projection  of  s  onto  lZ((Sn)'). 

Let  ||q||  >  0.  By  extending  to  the  present  case  (26),  one  finds  that,  as  h  — >  oo, 

Mh(x)  =  V{q)[l  +  0{h-1)]h-p^\\q\\  <Mh~l ||g||  (32) 

where  both  V(q)  and  p(q )  only  depend  on  (jf/||</||.  Further,  V(q)  and  M  are  positive  reals  and  p(q)  an 
integer  which  takes  values  in  p  ,  where  p  denotes  the  maximum  dimension  of  the  Jordan-blocks  associated 
to  the  Q-eigenvalues.  ■ 


In  the  sequel,  our  attention  will  be  focused  on  the  family  of  state-feedback  gain-matrices 

F:={Fh}™=h  ,  h>u  (33) 

along  with  the  system  (6)  under  a  time-varying  state-feedback  control  u(t)  =  Fhp\x{t ),  Fh^  G  F, 


x{t  +  1)  = 


$ 


h(t) 


:=  T  +  GF, 


(34) 


h(t) 
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4  Exponential  stability  under  arbitrary  admissible  switching 

A  control  horizon  sequence  {h(t)}t£z+,  h(t )  £  Z+,  h{t)  >  v,  is  called  admissible  if 

h{t  +  1)  >  h(t)  —  1  (35) 

Let  S  denote  the  set  of  all  such  admissible  sequences,  and  (£,<S)  the  system  (34)  under  an  arbitrary 
admissible  switching  sequence  in  S.  The  statements  that  close  the  present  section  are  proved  in  (Mosca, 
2005). 

Lemma  1  Along  all  possible  trajectories  of(Y,,S),  one  has 

\\x{t  +  h{t))\\h{t+h{t))  <  Mt)\\h{t)  (36) 

for  all  x(t)  such  that  ||s(i)||^m  >  0. 


The  theorem  that  follows  is  the  main  result  of  this  section  and  fundamental  for  the  subsequent  develop¬ 
ments  of  this  paper. 


Theorem  1  Consider  the  control  system  (£,<S)  composed  of  the  LTI  reachable  plant  (6)  under  a  time- 
varying  state-feedback  control  u(t )  =  Fh^x{t)  realized  by  arbitrary  admissible  switching  sequences  in  S. 
Then,  provided  that  h>  n  and  {h(t)}t£z+  be  bounded,  viz. 


n  <  h<  h(t )  <  h 

(37) 

with  h  <  oo,  (£,<S)  is  exponentially  stable 

\\x(t)  |  <  MXt  x(0) 

(38) 

with  M  S  (0,  oo),  and  decaying  rate  A  depending  on  h  and  h. 


Remark  3  -  It  is  to  be  pointed  out  that  (37)  encompasses  the  case  of  a  fixed  regulation  horizon  h, 
h  =  h  =  h  >  n,  for  which,  to  the  best  of  the  author’s  knowledge,  exponential  stability  of  (34)  for  an 
arbitrary  h  has  remained  so  far  an  open  problem  (Mosca,  1995).  ■ 

Remark  4  -  The  admissibility  condition  (35)  turns  out  to  be  not  only  sufficient  for  the  stability  property 
stated  in  Theorem  1,  but,  in  a  wide  sense,  also  necessary,  if  no  other  condition  is  imposed.  In  fact,  there 
are  cases  wherein,  if  (35)  is  violated,  stability  is  lost  even  if  h  remains  finite.  There  are,  of  course,  ways 
other  than  (35)  to  ensure  stability  under  arbitrary  state-gain  switching,  e.g.,  the  system  state  can  be 
suitably  processed  at  any  switching  time  as  indicated  in  (Hespanha  &  Morse,  2002).  ■ 
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5  Hysteresis  switching  regulation 


Noiseless  case  -  Consider  the  noiseless  variant  of  (1) 

x(t  +  1)  =  4>x(f)  +  Ga(u(t)) 


(39) 


along  with  (2)  and  (7).  Given  the  system  state  x,  Uh(x)  will  denote,  as  in  Sect.  3,  the  minimum  energy 
input  sequence  in  U^fx).  Consider  the  following  horizon  switching  logic  (h  >  n) 


h{t) 


h(t),  if  M-h{t){x{t))  <  1 
h(t),  otherwise. 


(40) 


h(t)  :=  max{4,  h(t  —  1)  —  1} 
h(t)  :  =  min{/i  e  Z+  :  h  >  h(t  —  1);  Mh(x(t))  <  1} 

with  t  =  1,2,...,  h(— 1)  =  h,  and  M/j( x)  as  in  (23).  Then,  the  following  theorem,  whose  proof  can  be 
found  in  (Mosca,  2005),  follows. 


Theorem  2  Consider  the  noiseless  input- saturated  system  (39),  (2)  and  (7)  with  u(t)  =  Fhnyx(t)  subject 
to  the  feedback  switching  logic  (40).  Assume  that  the  initial  system  state  x(0)  at  time  0  be  such  that 
h(0)  exist  finite.  Then,  logic  (40)  yields  the  admissible  switching  sequence  h(t  +  1)  =  h(t)  —  1,  t  = 
0, . . . ,  h(0)  —  h,  h{t)  =  h,  Vi  >  h( 0)  —  h,  the  resulting  switched  system  x(t  +  1)  =  ^hn^x(t)  satisfies 
the  input  saturation  constraints,  is  asymptotically  stable,  and  ||x(t)||  decays  to  zero  exponentially  fast.  In 
particular,  if  S  is  ANCBI,  the  resulting  switched  system  is  globally  asymptotically  stable  and  semi-globally 
exponentially  stable. 


Noisy  case  -  It  is  known  (Hou  et  al.,  1998;  Saberi  et  al.,  2000)  that  (27)  with  £<j(-)  =  0  has  the  most 
general  structure  of  an  input-saturated  LTI  system  for  which  it  makes  sense  to  consider  stability  and 
boundedness  under  arbitrary  loo-disturbances.  Under  such  circumstances,  the  argument  that  follows  is 
used  to  render  plausible  the  conjecture  that  h(t),  chosen  by  (a  suitably  modified  version  of)  (40),  cannot 
get  unbounded.  If  this  is  the  case,  according  to  Theorem  1,  the  contribution  of  any  past  noisy  sample 
vanishes  exponentially  fast,  and  hence  x(t)  stays  bounded. 

Consider  first  that,  by  the  switching  logic  (40),  h( 0)  <  oo  for  any  x(0)  S  M™,  and  {/i(i)}£f.0  Next, 

assume,  by  contradiction,  that  h(-)  grows  unbounded.  This  implies  that  ||x(-)||  does  the  same.  As  ||s(-)|| 
is  bounded  because  S  is  stable  and  a{u)  and  are  both  bounded,  there  are  times  t  large  enough  at  which 
llx(^)l|2  =  lls(^)l|2  +  ||g(t)||2  ~  ||g(t)||2.  Under  these  circumstances,  h(t)  is  essentially  chosen  according 
to  the  restricted  noiseless  system  with  state  q(t).  Hence,  according  to  Theorem  2,  at  subsequent  times, 
h(t  +  k)  =  h{t)  —  k  till  || q(t  +  k)  ||  is  reduced  so  much  that  the  effect  of  || s(t  +  A:)||  becomes  significant 
again  for  the  selection  of  subsequent  horizons.  Thereafter,  h(-)  may  start  to  increase  till  the  condition 
||x(f)||2  ~  ||g(t)||2  is  possibly  restored.  Consequently,  the  regulation  horizon  begins  again  to  decrease  by 
one  at  each  subsequent  time.  In  words,  a  “ horizon  resetting,'>  mechanism  is  inherently  enforced.  The 
conjecture  is  that  such  a  mechanism  prevents  h(-),  and  hence  the  plant  state,  from  becoming  unbounded. 
However,  in  order  to  prove  that  the  horizon  resetting  property  holds,  it  is  required  to  replace  the  switching 
logic  (40)  with  its  variant  (41)  equipped  with  a  “ hysteresis ”  facility.  The  latter  consists  of  the  hysteresis 
constant  e  in  (41)  which  ensures  that  the  values  taken  on  by  the  control  input  components  stay  strictly 
inside  the  saturation  levels  whenever  the  control  horizon  is  greater  than  h  and  required  not  to  decrease. 
Consequently,  e  >  0  in  (41)  makes  h  “easier”  to  decrease  than  to  increase  or  stay  constant. 
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Theorem  3  (Mosca,  2005)  Consider  a  generic  input- saturated  noisy  ANCBI  system  algebraically  equiv¬ 
alent  to  (27)  with  £q(-)  =  0.  Let  the  control  u(t )  =  Fhn\x(t )  with  Fh  as  in  (8)  and  h(t )  chosen  according 
to  the  following  hysteresis  switching  logic  (h>n) 


h(t) 


h(t),  if  M~Ht){x{t))  <  1 
h{t),  otherwise. 


(41) 


h(t)  :=  max{4,  h{t  —  1)  —  1} 

h(t)  :=  min{/i  G  Z+  :  h  >  h(t  —  1);  Mh(x(t ))  <  1  —  e} 

where:  t  =  1,  •  •  •;  e  G  (0,  1);  h(0)  =  h( 0)  with  h(— 1)  =  h;  and  Mh(z)  as  in  (23).  Then,  the  resulting 
closed-loop  hysteresis  switched  system  is  bounded-noise  bounded- state  loo-stable  irrespective  of  both  the 
initial  state  x(0)  G  Mn  and  the  magnitude  of  £<,(•). 


Proof.  See  the  Appendix. 


Theorem  3  ensures  loo-stability  irrespective  of  the  magnitude  of  the  state-disturbance  provided  that  the 
latter  does  not  affect  directly  the  unstable  modes  of  the  ANCBI  plant.  In  the  case  of  a  generic  state- 
disturbance  which  can  affect  directly  both  the  stable  and  the  unstable  modes  of  an  ANCBI  plant,  next 
Theorem  4  shows  that  the  ’’horizon  resetting”  mechanism,  and  hence  loo-stability,  holds  true  provided 
that  the  magnitude  of  the  state-disturbance,  directly  injected  on  the  critically  unstable  modes  be  smaller 
than  a  bound  dependent  on  the  hysteresis  constant  e. 

Theorem  4  (Mosca,  2005)  Consider  a  generic  input-saturated  noisy  ANCBI  system  agebraically  equiv¬ 
alent  to  (27)  with  £9(-)  /  0  and  u  =  =  Ui,  Vi  G  m.  Let  the  control  u(t )  =  Fh(t)x{t)  with  Fh  as  in 

(8)  and  h(t )  chosen  according  to  the  hysteresis  switching  logic  (41).  Then,  the  resulting  closed-loop  hys¬ 
teresis  switched  system  is  bounded-noise  bounded-state  loo-stable  irrespective  of  the  initial  state  x(0)  G  Mn 
provided  that,  Vt  G  Z+, 

UM\  <  e(l-e)/M  (42) 

where  e  is  the  hysteresis  constant  in  (41)  and  M  is  as  in  (26)  . 

Proof.  See  the  Appendix.  ■ 


Remark  5  -  The  switching  logic  (41)  ((40))  can  be  simplified  with  no  change  in  the  conclusions  of 
Theorem  3  and  Theorem  4  (Theorem  2)  by  trading-off  performance  vs.  memory  and/or  computational 
complexity,  e.g.,  by  replacing  h(t)  in  (41)  [(40)]  with 

h(t )  :=  rninjli  G  iU  {h(t  —  1)}  :  h  >  h(t  —  1);  Mh{x{t))  <  1  —  e}  (43) 

where  H  is  an  assigned  subset  of  Z+,  e.g.,  H  =  {h,h  +  A h,h  +  2A h,-  ■  •},  with  Ah  a  positive  integer. 
With  such  a  choice,  h(-)  can  only  increase  by  a  multiple  of  Ah  at  a  time.  On  the  other  hand,  according 
to  (35),  h(-)  can  only  decrease  by  one  unit  at  a  time.  In  this  connection,  one  could  avoid  to  hold  in  the 
computer  memory  all  feedback-gain  matrices  Fh,  Vh  G  {h,  ■  ■  ■,  h},  by  exploiting  the  fact  that,  as  can  be 
shown,  Uh-i(k\$h%)  =  Uh(k  +  l|x)  or  equivalently 

Fh.^k)  =  Fh(k  +  l)^1  (44) 

k  =  0, 1,  •  •  •  ,  h  —  2,  provided  that  &h  be  nonsingular.  ■ 
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Partial  state  information  -  If  in  the  hysteresis  switching  logic  (41)  the  true  plant  state  x(t)  is  replaced 
by  the  vector  x(t)  +  [£7(f)  O']7  where  £(t)  G  Mns  is  a  bounded  sensor-noise  acting  on  the  stable  component 
of  the  x-state,  it  is  immediate  to  see  that  the  conclusions  of  Theorem  3  hold  true.  This  implies  that  the 
statement  of  Theorem  3  can  be  extended  to  an  ANCBI  plant  algebraically  equivalent  to  (27)  under  the 
hysteresis  switching  logic  (41)  based  on  a  partial  state  information,  x(t)  =  [s'(i)  q'(t)]'  where  s(t)  is  a 
filtered-estimate  of  s(t)  based  on  observations  y(t)  =  [y(t)  q'(t)]'  with  7 (t)  =  Hs(t )  +  n(t)  E  Mp  with  n(-) 
a  bounded  noise.  More  generally,  if  a  bounded  sensor-noise  acts  also  on  the  unstable  component  of  the 
2-state,  Theorem  4  can  be  used  so  as  to  estimate  the  noise  magnitude  below  which  boundedness  can  be 
guaranteed. 

6  An  Example 

The  reader  is  referred  to  (Bacconi,  Casavola  &  Mosca,  2003)  for  more  details  on  the  example  reported  here¬ 
after  concerning  a  multiple-input  satellite  formation  problem.  Formation  flying  of  satellites  is  currently  an 
active  area  of  research.  The  simplest  problem  in  formation  flying  control  is  the  formation  keeping  in  which 
the  formation  is  composed  by  two  spacecrafts.  In  such  a  case,  the  so  called  “master/slave”  approach  can 
be  used  to  maintain  a  desired  position  of  the  slave  satellite  with  respect  to  the  master  satellite.  Because 
the  input  forces  are  subject  to  saturation,  MBPC  strategies  can  be  considered  particularly  effective  to 
that  end.  A  standard  approach  to  the  problem  of  formation  flying  control  is  to  linearize  the  dynamics  of 
the  satellites  around  some  reference  orbit.  In  the  case  that  the  reference  orbit  is  approximately  circular 
and  for  small  distances  (less  than  1  km)  between  the  two  vehicles,  the  orbit  and  position  of  a  pair  of 
nano-satellites  can  be  schematized  as  in  Fig.  6.  From  the  related  linearized  equations  of  motion  we  get 


Figure  1:  Position  of  the  slave  satellite  with  respect  to  the  master. 


the  following  state  space  equations  of  the  relative  motion  of  the  slave  spacecraft: 

x  =  4>cx  +  Gc(u  +  d) 


where 


x  =  [xp  xp  yp  i/p}' 
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Here:  xp,  yv  are  the  relative-coordinates  of  the  slave  satellite;  u  =  [ux  uy  ,  ]7  are  the  actuator  forces  in  the 
xp  and  yp  positive  axis  directions;  d  =  [dXp  dVp]'  are  the  relative  disturbances.  Further,  co  =  0.001086  rad/s 
is  the  assumed  angular  velocity  of  the  master  satellite  in  the  inertial  frame.  With  a  sampling  period  T, 
the  ZOH  sampled  dynamical  model  of  the  slave  satellite  takes  the  form 

x(t  +  1)  =  <hx(t)  +  Gu(t)  +  £(t)  (45) 

In  (45),  £(t)  represents  relative  disturbance  forces  accumulated  during  a  sampling  period  of  duration  T. 
At  any  sample  time,  a  measure  x(t)  of  the  state  is  available:  x(t)  =  x(t)  +  x(t)  with  x(t)  the  measurement 
error.  The  objective  of  the  control  is  to  drive  the  slave  vehicle  to  the  desired  relative  position  (xp,yp). 
Subsequently,  the  relative  position  has  to  be  maintained  in  spite  of  the  presence  of  disturbances  and  state 
measurement  errors,  with  an  accuracy  of  at  least  1  cm.  Manoeuvres  are  made  thanks  to  a  combination 
of  small  jet  actuators  able  to  produce  any  force  components  within  a  given  square  rectangular  region 
U  =  [— um,um]  x  [— um,um\.  Hence,  equation  (45)  has  to  be  replaced  by 

x(t  +  1)  =  <&x(i)  +  Ga(u(t))  +  £(t) 

We  apply  the  control  procedure  of  this  paper  to  the  problem  of  formation  flying  control  of  a  pair  of  nano¬ 
satellites  of  masses  equal  to  10  kg.  We  consider  a  sampling  period  T  =  0.3  s.  Note  that  the  discrete-time 
system  (45)  is  ANCBI.  The  disturbance  set  is  assumed  to  be  spherical  with  radius  £o  =  0.0025.  The 
objective  of  the  control  is  to  keep  the  slave  craft  in  a  ball  set  centered  on  (xp,yp)  =  (100,0)  with  radius 
0.01  m.  The  input  set  U  is  assumed  to  be  an  axis-aligned  rectangle  with  opposite  corners  at  [—0.15,  —0.15] , 
[0.15,  0.15]. 

The  following  figures  show  simulation  results  for  an  initial  condition  x(0)  =  [57  3  17  l]r.  Stability  and 
positioning  accuracy  is  good  in  spite  of  the  large  (<  0.002)  random  errors  in  the  state  vector  estimates. 
We  observe  that  the  desired  position  of  the  slave  satellite  is  obtained  in  less  than  150  time  steps,  that  is 
about  45  s  (Fig.  3(a)).  Thereafter,  the  formation  keeping  is  achieved,  with  an  accuracy  of  less  than  1  cm. 
Fulfillment  of  the  input  saturation  constraints  is  highlighted  by  Fig.  4(a).  The  results  are  comparable  or 
better  than  others  obtained  by  alternative  techniques.  The  evolution  of  the  regulation  horizon  is  depicted 
in  Fig.  4(b).  Notice  that,  because  of  the  initial  state  and  the  input  saturation  constraints,  the  initial 
value  of  the  required  regulation  horizon  is  large:  about  150  time  steps.  It  guarantees  that  all  the  predicted 
samples  of  the  input  sequence  satisfy  the  constraints.  Afterwards,  the  value  of  h(t)  decreases  by  one  at 
each  sampling  time  and  reaches  the  dimension  of  the  system. 
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Figure  2: 


(a)  Evolution  of  the  relative  position  of  the  slave  (b)  Trajectory  of  the  slave  satellite  with  respect 
satellite  with  respect  to  the  master.  to  the  master. 


Figure  3: 
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(a)  Input  forces  to  the  slave  spacecraft. 
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7  Conclusions 

The  paper  provides,  relatively  to  alternative  approaches,  a  computationally  affordable  solution  to  the  reg¬ 
ulation  problem  of  discrete-time  input-saturated  LTI  systems  subject  to  persistent  bounded  disturbances 
of  unknown  arbitrary  magnitude.  The  proposed  solution  enjoys  the  following  features:  It  consists  of  a 
supervisory  switching  control  logic  whereby  a  feedback-gain,  selected  at  any  time  from  a  family  of  pre¬ 
designed  candidate  feedback-gains,  is  switched-on  in  feedback  to  the  plant  according  to  the  information, 
either  complete  or  partial,  on  the  current  plant  state;  No  disturbance  upper-bound  need  to  be  known; 
The  controller  selection  is  made  in  accordance  with  a  predictive  control  philosophy,  and  each  candidate 
feedback-gain  is  tuned  on  to  a  different  horizon  in  a  receding-horizon  control  sense;  The  supervisory 
switching  logic  is  flexible  enough  so  as  to  enable  the  designer  to  simplify  the  scheme  by  trading  off  perfor¬ 
mance  vs.  memory  and/or  computational  complexity  while  retaining  the  required  stability  properties.  It 
is  proved  that  the  adopted  switching  logic  ensures  global  asymptotic  and  senri-global  exponential  stability 
in  the  ideal  noiseless  case,  and  finite  -induced  gain  to  the  disturbance-to-state  map,  whenever  the 
structure  of  the  disturbed  plant  can  make  such  properties  conceptually  achievable,  viz. ,  the  disturbance 
which  enters  an  ANCBI  system  does  not  act  directly  on  the  critically  unstable  modes,  the  latter  being 
only  indirectly  affected  by  the  disturbance  via  the  feedback  controls.  If  this  is  not  the  case,  a  hysteresis 
facility  in  the  switching  logic  becomes  essential  for  handling  disturbances  of  small  enough  magnitude 
which  directly  affect  the  critically  unstable  modes  of  the  plant. 

For  unstable  (non- ANCBI)  noisy  systems  with  input  saturations,  finite-gain  stabilization  is  a  rather  in¬ 
tricate  matter.  In  such  a  case,  it  is  known  that  no  global  or  semi-global  finite-gain  stabilization  can  be 
achieved.  In  particular,  suitability  of  our  approach  for  that  case  appears  yet  unexplored  as  in  the  present 
paper  our  methods  of  proof  mainly  rely  on  property  (22)  whose  validity  is  limited  to  ANCBI  systems. 
The  approach  adopted  in  this  paper  lies  at  the  crossroad  of  different  control  methods:  constrained  MBPC 
whereby  an  optimization  problem  is  numerically  solved  in  real-time  within  each  sampling  interval;  con¬ 
straint  management  via  a  reference-governor  unit  (Gilbert  &  Kolmanovsky,  2002)  whereby  a  fixed  feedback 
compensator  is  used,  and,  hence,  the  on-line  optimization  is  restricted  to  only  re-modulate  the  shapes 
of  the  set-point  changes;  and  switching  supervisory  control.  In  this  paper  robustness  against  persistent 
bounded  disturbances  of  unknown  arbitrary  magnitude  is  acquired,  in  contrast  with  standard  MBPC,  at 
the  cost  of  giving  up  on-line  optimization  (and,  hence,  full  performance),  while  circumventing  the  bot¬ 
tleneck,  present  in  reference  governor  methods,  of  a  fixed  preselected  compensator,  by  switching  among 
precomputed  candidate  linear  feedback-gains. 
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